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ABSTRACT 

We report the detection of [Ne n] emission at 12.81 fim in four out of the six opticaUy thick 
dust disks observed as part of the FEPS Spitzer Legacy program. In addition, we detect a 
H i(7-6) emission hne at 12.37/im from the source RX J1852. 3-3700. Detections of [Ne n] hnes 
are favored by low mid-infrared excess emission. Both steUar X-rays and extreme UV (EUV) 
photons can sufficiently ionize the disk surface to reproduce the observed line fluxes, suggesting 
that emission from Ne"*" originates in the hot disk atmosphere. On the other hand, the H i(7- 
6) line is not associated with the gas in the disk surface and magnetospheric accretion flows 
can account for at most ~30% of the observed flux. We conclude that accretion shock regions 
and/or the stellar corona could contribute to most of the H i(7-6) emission. Finally, we discuss 
the observations necessary to identify whether stellar X-rays or EUV photons are the dominant 
ionization mechanism for Ne atoms. Because the observed [Ne n] emission probes very small 
amounts of gas in the disk surface (~ 10~^Mj) we suggest using this gas line to determine the 
presence or absence of gas in more evolved circumstellar disks. 

Subject headings: line: identification - circumstellar matter - planetary systems: protoplanetary disks 
- infrared: stars - stars: RX Jllll.7-7620, PDS 66, HD 143006, [PZ99] J161411. 0-230536, RX J1842.9- 
3532, RX J1852.3-3700 



1. Introduction 

Young stars are often surrounded by gas and 
dust disks that may succeed in forming planets. 
The properties and evolution of their dust and gas 
components are key to understanding planet for- 
mation and the diversity of extrasolar planetary 
systems. 

Circumstellar dust has been extensively stud- 
ied in young and old disks since the IRAS mis- 
sion. Grain growth has been identified in disks 



around star s with masses ranging from a few so- 



lar masses (Bouwman et al. 2001) down to the 



brown dwarf regime (|Apai et al.ll2005h suggesting 
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that the first steps of planet formation are ubiq- 
uitous. Detailed mineralogy reveal chemical pro- 
cessing similar to those that occurre d in the early 
solar system (see iNatta et al.l 12007 for a review 
on du st i n protoplanetary dis ks; Alexander et al.l 
20071 and lWooden et al.ll2007l for dust processing 
in the solar nebula and in disks). On the other 
hand, the properties and evolution of circumstel- 
lar gas are less well characterized. 

Filling this gap is one of the goals of the Forma- 
tion and Evolution of Planetary Sy stems (FEPS) 
Spitzer legacy program (jMever et al.. 200 61 . In two 
previous contributions we derived gas mass upper 
limits for a sample of 16 sun- like stars surrounded 
by optically thin dust disks and discussed impli- 
cations for the formation of terrestrial, giant, and 
icy p lanets ( HoUenbach et al.ll2005t Pascucci et al 



20061) . 
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Here, we present high-resolution Spitzer spec- 
tra for the six FEPS targets with excess emis- 
sion beginning at or shortward of the 8 
IRAC band. Their IRAC colors are consistent 
with those of ac c reting classical T Tauri stars 



( Silverstone et al. 



20061) suggesting that these 
stars are surrounded by optically thick dust disks. 
We report the detection of the H I (7-6) line at 
12.37 /im in one system and of the [Ne li] line 
at 12.81 /im in four systems (Sect. 13. ip . Be- 
cause Ne atoms have a large ionization potential 
(21.6eV), the detection of [Ne ii] hnes is of par- 
ticular interest to assess the role of stellar X-ray 
and EUV {hv > 13.6 eV) photons on the disk 
chemistry and to explore the conditions for disk 
photoevaporation dGlassgold. Naiita fc Igeall2007 



iGorti fc HollenbachI 120071 ). We discuss in detail 
predictions from the proposed X-ray and EUV 
models and future observations that will be able 
to identify the dominant ionization mechanism 
(Sect.lHand Sect. [5]). Since both models find that 
the [Ne il] line probes small amounts of gas on 
the surface of circumstellar disks, we also suggest 
using this tracer to place stringent constraints on 
the presence or absence of gas in more evolved 
disks. 

2. Observations and Data Reduction 

The general properties of the 328 stars in th e 
FEPS sample are described in lMever et al. I (|2006h . 
We summarize in Table [T] the main properties 
of the six FEPS sources surrounded by optically 
thick dust disks. In Cols. 4, 5, and 6 we give the 
star spectral types (SpTy), distances (d), and ages 
(Age). References for these quantities are pro- 
vided in Col. 7. The stellar effective temperatures 
(Tefi), visual extinctions (^v)i Sind bolometric lu- 
minosities (i*) are listed in Cols. 8, 9, and 10. 
The last column summarizes the source X-ray lu- 
minosities (Lx)- 

To compute the X-ray luminosities we used 
the ROSAT PSPC All-Sky Survey coun t -rates 



and HRl hardn e ss ra t ios ( Alcala et al.l 19971: 



Voges et al.l Il999l l2000[ ) following [Fleming et al.l 

and adopting the distances in Table [1] 
These X-ray luminosities are representative for 
the energy band 0. 1-2.4 keV (or 120-5 A). The er- 
rors in Lx include the uncertainties in the count- 
rates, HRl, and distances. Note however that the 



intrinsic variability of log(Lx) due to stellar ac- 
tivity is larger than the quoted uncertainty and 
amounts to at lea st a few tenths of a dex (e.g. 
Marino et al.ll2003h . Our X-ray luminosities agree 



with values fr om the literature ('Alcala et al. 1997 
for source 1; iMamaiek et al.l l2002 f or source 2; 



[Sciortino et al.lll998l for source 3; iNeuhauser et al 
2000l for sources 5 and 6). Because source 4 ap- 



pears extended in the ROSAT PSPC image, we 
also searched for an independent measurement of 
its X-ray flux. Its count-rates from the XMM- 
Newtoifl MOSl and M0S2 cameras in the 0.2- 
2keV bandpass convert to log(Lx)=30.8 erg/s. 
Since an increase of 0.26 dex in luminosity could 
be simply due to stellar activity, we prefer to adopt 
the luminosity from ROSAT for consistency with 
the other sources. 

Among the sources listed in Table [l] HD 143006 
is the only one that was included in the FEPS pro- 
gram as part of the gas detection experim ent be- 



ing a known dust disk from IRAS (Svlvester et al 
19961 ). The other five dust disks have been re- 



cently identified by our group from IRAC / Spitzer 
photometry ( Silverstone et al.l 2006). The spec- 
tral energy distributions of these five systems, in- 
cluding photometry from IRAC and MIPS as well 
as IRS low-reso l ution spectra, are presented in 
Silverstone et al. ( 20061 ). Hillenbrand et al. in 
prep, show that single temperature blackbody fits 
to the 33 and 70 /im data cannot account for the 
excess emission at shorter wavelengths in any of 
the six sources in Table [U indicating the pres- 
ence of warm inner disk material. These optically 
thick dust disks are also the only six FEPS sources 
exhibiting dust fe atures in the Spitzer/I RS low- 
resolution spectra. iBouwman et al. (2007) present 
a detailed analysis of their mineralogy and find 
that RX J1842. 9-3532 is surrounded by an almost 
primordial disk (flared geometry and 1 /im~sized 
grains) while [PZ99] J161411. 0-230536 has the 
most processed disk (close to flat geometry and 
large 5/im grains). 

In Sect. 12.11 we describe the observations and 
data reduction of the high-resolution infrared 
spectra. We complemented the Spitzer data with 
optical spectra (Sect. 12. 2p that are used to esti- 



mate mass accretion rates (Sect. 



iXMM-Newton Serendi pitous Sour ce Catalogue, IXMM, 
|http : / / xm mssc- www , star . le . ac . uk/ 1 
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Table 1 

Main properties of the targeted optically thick dust disks. 



ID# 


Source 


2MASS J" 


SpTy 


d 


Ago 


Rof. 


Totf 


A, 


log(L*) 


log(Lx)'' 










[pc] 


[Myr] 




[K] 


[mag] 


[Lq] 


[erg/s] 


1 


RX Jllll. 7-7620 


11114632-7620092 


Kl 


163±10 


5 


1,2,3 


4621 


1.30 


0.27 


30.56±0.16 


2 


PDS 66 


13220753-6938121 


Kl 


86±7 


17 


4 


5228 


1.22 


0.10 


30.20±0.09 


3 


HD 143006 


15583692-2257153 


G6/8 


145±10 


5 


5,6,7 


5884 


1.63 


0.39 


30.14±0.17 


4 


[PZ99] J161411. 0-230536 


16141107-2305362 


KO 


145±10 


5 


8,6,7 


4963 


1.48 


0.50 


30.54±0.10 


5 


RX J1842. 9-3532 


18425797-3532427 


K2 


140±10 


4 


9,10,3 


4995 


1.03 


0.05 


30.37±0.18 


6 


RX J1852. 3-3700 


18521730-3700119 


K3 


140±10 


4 


9,10,3 


4759 


0.92 


-0.17 


30.57±0.13 



^The 2MASS source name includes the J2000 sexagesimal, equatorial position in the form: lihmmssss+ddmmsss JCutri et al.|[2003l) . 
^X-ray luminosities in the ROSAT 0.1— 2.4kcV energy band, see Sect.[2]for details. 



Note. — Spectral types are from optical spectroscopy with accuracy of 1 subtype. The stellar colors B— V and V-K are used in conjunction 
with the spectral type to estimate effective temperatures (Tgff). Visual extinctions (Av) are computed from the spectral type and stellar 
colors. We refer to Carpenter et al. 2007 in preparation for details. Stellar luminosities are computed from the best fit Kurucz stellar 
models to optical and near-infrared observations of the stellar photosphere (Carpenter et al. 2007 in prep.) and the distances in this table. 

Reference s. — (1 ) [Alcala et al. 1995; (2) Luhman 200 ^; (3) Hillenbra nd e t al. 2007 in prep.; (4) [Mamaick ct al. 2002; (5) 
[Hou k fc Smith-Moore||l98&; (6) Ac Zeeuw et al.. 1999.: (7) Preibisch et al.|f200l (8) IPreibisch et al.|fT998l ; (9) [Neuhauser et al...2000. ; (10) 
iNeuhauser & Forbrichll20Q7l 



2.1. High resolution IRS Spectra 

Spitzer/IRS high-resolution spectra for the 6 
FEPS optically thick dust disks were obtained be- 
tween August 2004 and September 2005. Observa- 
tions were done in the Fixed Cluster-Offsets mode 
with two nod positions on-source (located at 1/3 
and 2 /3 of the slit length) and two additional sky 
measurements (1' east of the nodi and nod2 posi- 
tions) acquired just after the on-source exposures. 
We used these sky exposures to remove the in- 
frared background an d reduce the number of rogue 
pixels as described in IPascucci etaP (|2006l ). The 
PCRS or the IRS Peak-up options were used to 
place and hold the targets in the spectrograph slit 
with positional uncertainties always better than 
1" (1 sigma radial). Exposure times were cho- 
sen to detect a 5% line-to-continuum ratio with 
a signal-to-noise of 5 (see Table ^ . In the fol- 
lowing we focus on the reduction and analysis of 
the SH module (9.9-19.6 ^m) where we detected 
gas lines in four out of six targets. No gas lines 
are detected in the wavelength range 18.7-37.2 /im 
that is covered by the LH module. The SH mod- 
ule is a cross-dispersed echelle spectrograph with 
a resolving power of 650 in the spectral range 
from 9.9-19.6/im, corresponding to a spectral res- 
olution of 0.015 /im around 10 ^m. The detector 
has a plate scale of 2.3"/pixel and the slit aperture 
has a size of 2 x 5 pixels, thus covering a region of 
~ 1600 X 700 AU around a star at 140 pc. 



Da. ta reduction was carried out as in lPascucci et al 
( 20061 ) starting from the Spitzer Science Center 
(SSC) S13.2.0 droop products. We fixed pixels 
marked bad in the bmask files with flag value 
equal to 2^ or larger, thus including anomalous 
pixels due to cosmic-ray saturation early in the 
integration, or preflagged as unresponsive. We 
also inspected visually all the SH exposures to 
catch additional rogue pixels and found less than 
5 per frame. These bad and rogue pixels were cor- 
re cted using the SSC irs dean package as explained 



Pascucci et all (|2006f) 



We flux calibrated the extracted spectra using 
nine independent observations (over four differ- 
ent Spitzer campaigns, from C21 to C24) of the 
bright standard star ^ Dra and the MARCS stel- 
lar atmosphere model degraded to the spectro- 
graph's resolution and sampling. The dispersion 



in the mean fluxes of these calibrators is always 
less than 1.5% between 10 and 20 /im. We veri- 
fied the absolute flux calibration by comparing the 
high-res olution spectra to the low -resolution spec- 
tra from Bouwman et al. ( 2007t ) at a wavelength 
in between the [Ne ii] and the H i(7-6) lines. We 
found deviations less than ~10% at 12.6±0.1/im, 
that are within the absolute photometric accuracy 
estimated by the SSC, for all objects except for 
RX Jllll. 7-7620. For this source the SH spec- 
trum has 30% higher flux but maintains the same 
shape as the low-resolution spectrum. We verified 
that this higher flux is simply due to the non-zero 
background emission in the vicinity of the source 
(the sky position we acquired for the SH spectrum 
is too far away from the source to be representative 
of the nearby emission). We have also tested that 
the low-resolution spectra agree with the IRAC 
photometry at 8 /xm and the MIPS photometry at 
Thus, to improve our flux estimates and 
upper limits we scaled the high-resolution spectra 
to the flux of the low-resolution spectra at 12.6 /im. 

2.2. Optical Spectra 

Optical spectra were obtained in several observ- 
ing runs and with different telescopes. The spectra 
of HD 143006, RX J1842.9-3532, and RX J1852.3- 
3700 were acquired in July 2001 and June 2003 
with the Palomar 60-inch teles cope and facilit y 
spectrograph in its echelle mode (|McCarthvlll988h . 
A 1'.'43 X 7'.'36 slit was used in combination with a 
(2-pixel) resolving power of ~19,000. Details on 
the observational strategy an d data reduc t ion ar e 
provided in Sects. 2 and 3 of White et al. ( 20061 ). 
White et al] (|2006[ ) also derive several stellar prop- 
erties from the Palomar 60-inch spectra for these 
and other FEPS sources, including radial and ro- 
tational velocities, Li i A 6708 and Ha equivalent 
widths, chromospheric activity index R'^xj ^'^^ 
temperature- and gravity-sensitive line ratios. 

The spectrum of [PZ99] J16141 1.0-230536 was 
acquired during commissioning of the East Arm 
Echelle spectrograph on the Hale 200-inch tele- 
scope at Palomar. The on-source exposure time 
was 2,000 seconds. The spectrum covers most 
of the 4000-9350 A region with a resolution that 
varies from 16,000 to 27,000 within each order. We 
extracted the spectrum and corrected for scattered 



^http: / /ssc. spitzer. caltcch.edu/irs/calib/tcmpl/ 



^http: / /ssc. spitzer. caltech.edu/legacy/fepshistory.html 
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Table 2 

Log of the IRS short-wavelength high-resolution observations. 



Source 


AOR Key 


SH 


Peak-up 






time X ncyclcs 


mode 


RX Jllll. 7-7620 


5451776 


31.5x12 


PCRS 


PDS 66 


5451264 


6.3x4 


PCRS 


HD 143006 


9777152 


6.3x18 


PCRS 


[PZ99] J161411. 0-230536 


5453824 


31.5x4 


PCRS 


RX J1842. 9-3532 


5451521 


31.5x10 


IRS-red 


RX J1852. 3-3700 


5452033 


121.9x4 


PCRS 



light using custom rout ines in IDL, following the 
main steps described in White et al.l ([2006). 

The spectrum of RX Jill 1.7- 7620 was ac- 
quired in March 2003 with the MIKE echelle 
spectr ograph on the Magell an Clay 6.5-m tele- 



scope (jBernstein et al.M2002[ ). The star was ob- 



served in a 360~second exposure with the MIKE 
Red CCD in the standard setup with the 0.35" 
slit, and 2-pixel resolution of R~ 36,000 covering 
the wavelength range 4800-8940 A. The data were 
reduced using the MIKE Redux IDL package^. 

Finally, the spectrum of PDS 66 was acquired 
in April 2002 with the echelle spectrograph on the 
CTIO 4 meter Blanco telescope. We used the 31.6 
red long echelle grating covering the wavelength 
range between 3,000-10,000A and a 0'.'8x3'.'3 slit 
with a 2-pixel resolution of R~45,000. The ex- 
posure time on-source was 120 seconds. For 
the data reduction we used the IRAF pack- 
ages quadred I echelle that can treat multiamplifier 
echelle data. 

3. Results 



In Sect. 13.11 we present the gas line detections 
from the Spitzer spectra of the six FEPS optically 
thick dust disks. We also compute mass accretion 
rates from the Balmer emission profiles (Sect. 13. 2|) 
and use them in Sect. |4] to search for correlations 
with the flux of the detected infrared lines. 

3.1. Fluxes of Infrared Gas Lines 

We report the detection of [Ne ii] emission at 
12.81 /zm in four out of the six targets and the 
additional detection of H i(7-6) at 12.37 ^m in 
RX J1852.3-3700 (see Fig. [T]). For these detec- 



tions, we verified that the emission is centered at 
the expected rest wavelength, is spectrally unre- 
solvecO, is present in both nodi and nod2 posi- 
tions, and is absent in the sky positions. These 
checks guarantee that the emission originates at 
the source location or in its proximity, within 
2'.'3/pixel or ~320AU at a distance of 140 pc. 

To estimate the fiuxes of the detected lines, we 
fit the spectrum within ±0.25 fim of each line using 
a Levenberg-Marquardt algorithm and assuming 
a Gaussian for the line profile and a first-order 
polynomial for the continuum. The la errors on 
the line fiuxes are evaluated from the RMS dis- 
persion of the pixels in the spectrum minus the 
best fit model (see Table ^ . In the case of non- 
detections, we fit the same spectral range with a 
first-order polynomial and we provide in Table [3] 
the 3ct upper limits to the flux from the RMS in 
the baseline subtracted spectrum. Figs. [1] and [2] 
illustrate our best fits to the detected lines (red 
dashed lines) and the hypothetical 3(T upper lim- 
its (light blue dot-dashed lines). In addition to the 
H I and the [Ne ii] transitions, we provide upper 
limits for the [Ne ill] line at 15.55 /im for which we 
have predicted line emissions from the X-ray and 
EUV models. 

3.2. Mass Accretion Rates 

Observational evidence for magnetospheric ac- 
cretio n in classical T Tau ri stars (TTSs) is robust 
(e.g. Bouvier et al.l 2007 ). Evidence of this pro- 
cess are the UV/optical continuum excess, emit- 
ted by the accretion shock on the stellar surface, 
and permitted emission lines origi nating in the in 



falling magnetospheric gas (e.g. iMuzerolle et al 



*http:/ /web. mit.edu/ burles/www/MIKE 



®The resolving power of the SH module is R~650 corre- 
sponding to '~460km/s. 
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Fig. 1. — Expanded view of the wavelength regions around the H i(7-6) and [Ne ii] emission hnes. On top 
of the stellar and dust continuum we overplot the best Gaussian fits to the data (red dashed-lines) and the 
hypothetical 3a upper limits (light blue dot-dashed lines) reported in Table [31 In the case of PDS 66, we 
might have detected [Ne ii] emission at a level of ~2cr (Flux^^ 1.4 x 10~^^ergs~^cm~^). However, due to 
the faintness of the emission we cannot confirm its presence in both nod positions and therefore prefer to 
report a 3a upper limit in Table O 
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Fig. 2. — Expanded view of the wavelength regions around the [Ne ill] Hne at 15.55 /im. On top of the steUar 
and dust continuum we overplot the hypothetical 3ct upper limits (light blue dot-dashed lines) reported in 
Table [1 



Table 3 

Line fluxes and upper limits (3(t) from the IRS/Spitzer spectra. 



Source 


Flux(H i) 

[10^^^ crgs^^cm"^] 


Flux([Nc ii]) 


Flux([Ne III]) 
[10^^^ orgs"^cm"^] 


RX Jllll. 7-7620 


<5.6 


5.1±1.2 


<3.0 


PDS 66 


<26.0 


<20.0 


<21.1 


HD 143006 


<22.6 


<15.0 


<19.7 


[PZ99] J161411. 0-230536 


<7.2 


6.2±1.8 


<4.7 


RX J1842. 9-3532 


<6.1 


4.3±1.3 


<3.6 


RX J1852. 3-3700 


3.4±0.4 


7.2±0.4 


<1.5 



Note. — The wavelengths for the H i(7-6), [Ne n], and [Nc ill] Unes are 12.37, 12.81, and 15.55 /^m 
respectively. The Icr errors on the detections are calculated from the RMS of the observations minus 
model fit. 
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20011: iKurosawa et al.|[2006l ). The UV/optical con- 
tinuum excess emission is the most common and 
direct measure of mass accretion rates (M*). How- 
ever, in many instances, this excess is too weak to 
be measured, and other methods such as model- 
ing Balmer emission profiles are necessary. This is 
particularly true for objects with predominately 
low accretion rates, such as older 5-10 Myr TT Ss 
(e.g. iMuzerolle eraLlbOOd iLawson et al.ll2004 ) in 
the same age range as our sample. 

Many diagnostics in our spectra suggest low ac- 
cretion rates in comparison to values typical of 
younger TTSs such as in Taurus: the lack of opti- 
cal continuum veiling (continuum excess over pho- 
tospheric flux <0.1), weak or absent mass loss 
signatures such as [O i] A6300A emission, and 
the lack of broad emission components in the 
Na D doublet and Ca ii triplet lines. The Ha 
and H/3 emission profiles, shown in Figures [3] and 
m are in most cases qualitatively suggestive of 
weak accretion. One object, RX J1842. 9-3532, ex- 
hibits blueshifted absorption indicative of signifi- 
cant mass loss, and two objects exhibit redshifted 
absorption components in Ha ([PZ99] J161411.0- 
230536 and HD 143006). Given the spectral 
types of these stars, the lack of optical veiling im- 
plies an upper lim it of M^, < 10~^ Mq yr"-'^ (see 
Calvet et al1l2004[ ). 



We followed the procedures outlined in detail 



Muzerolle et al 



( 200lh to model the Balmer 
line profiles and thus estimate mass accretion 
rates for our targets. Although the paper by 



targi 

Muzerolle et al.l (|2001.) focused on low-mass stars 



with late K and M spectral types, ac cretion diag 



nosti c s such as UV exce ss and Pa/3 (jCalvet et al 



2004 . Natta et al. 20061 ) indicate that magneto- 



spheric model assumptions can be extended to 
stars with early K and G spectral types like our 
targets. We adopted stellar parameters based on 
the empirically-derived quantities in Table [TJ Gas 
temperatures were set roughly to 10,000 K fol- 
lowing the constraints derived in IMuzerolle et al 



(|200l[ ). Note that for the density regime appropri- 
ate for these particular objects, the gas is already 
nearly fully ionized and larger gas temperatures 
will not result in significantly different line emis- 
sion. The magnetosphere size/ width was set to 
a fiducial value (between 2.2-3 i?*), lacking any 
observable constraints. This is the largest source 
of uncertainty in constraining the accretion rate; 



however, a plausible range of values results in no 
more than a factor of 3-5 range in accretion rates 
that can reproduce the observed line profiles. We 
further included rotat ion, using the treatment of 



Muzerolle et al.l (|2001l ). since rotation rates typi- 



cal of our objects can have an observable effect on 
the line profiles near the line center. The adopted 
stellar equatorial velocities, based on the observed 
V sin{i) values, are given in Table [H Finally, the 
model inclination angle and mass accretion rate 
are varied to reproduce the wings, width, and peak 
of the observed lines. The best matches are shown 
in Figs. [3] and [4] and the inferred inclination angles 
(i) and mass accretion rates (M^,) are summarized 
in Table H Note that < 5 x 10~" M© yr'^ 
produce negligible Ha emission for our sources 
and thus represent lower limits for their accretion 
rates. Given the magnetospheric model param- 
eters in Table IH we also compute the predicted 
flux in the H i(7-6) transition (last column of the 
table). 

The models reasonably account for the observed 
profiles, with a few exceptions. The Ha emission 
line of RX J1852. 3-3700 exhibits a narrow, sym- 
metric core that the models cannot reproduce. It 
is possible that there is a chromospheric emission 
component superposed on top of the broader ac- 
cretion component. For this to occur, a significant 
amount of the stellar surface must not be covered 
by the accretion fiow (which would be consistent 
with the more pole-on orientation of the model). 
The models also cannot reproduce the blueshifted 
absorption components seen in both Ha and H/3 
profiles of RX J1842.9-3532. A treatment of the 
wind is necessary to account for such features, 
which is beyond the scope of this work. How- 
ever, we note that the models still account for the 
overall emission profile, suggesting that the wind 
produces negligible emission. 

The derived mass accretion rates are all lower 
than the average value for 1 Myr-old T Ta uri stars 
(|Gullbring et al.lll998l : ICalvet et alll2004l) by fac- 
tors larger than 10. Given the older ages of the 
stars in our sample, this may be consistent with 
a general decline in accretion rate with time as 
predicted by models of viscous disk evolutio n (e.g. 
Hartmann et al.lll998l:rMuzerolle et~aIl l2000V The 



predicted H i(7-6) emission for RX J1852.3-3700 
contributes to at most '--^30% of the observed 
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valuti, while predicted fluxes from the other 
sources are negligible or fall well below the So- 
upper limits we report in Table [H Thus, mag- 
netospheric accretion models argue against the 
H i(7~6) transition originating in accretion flows. 
Accretion shock regions, the stellar chromosphere, 
and the hot disk surface are other possible sources 
of emission. We discuss the contribution from the 
disk in the next Section and note that the possi- 
ble chromospheric emission component in the Ha 
profile of RX J1852. 3-3700 could also contribute 
to the H i(7-6) line. 

4. Disk Atmosphere and [Ne ii] Emission 

Two-thirds of our optically thick dust disks 
present a [Ne ii] emission line at 12.81 /im. Where 
does the [Ne ll] emission orig inate? Recently 
Glassgold. Naiita fc Igeal (2007) showed that stel- 
lar X-rays can partially ionize the gas in the disk 
atmosphere and produce detectable [Ne ii] lines. 
Alternatively, HoUenbach & Gorti in prep, sug- 
gest that extreme ultraviolet (EUV, hi/ > 13.6 eV) 
photons from the central star ionize the upper 
layer of circumstcUar disks and create a kind of 
coronal H ii region producing detectable [Ne ii] 
emission. Although the two models rely on dif- 
ferent ionization mechanisms, they both predict 
that [Ne ii] emission originates from gas in a hot 
surface layer of circumstellar diskfl We can thus 
expect to find some correlations between the ob- 
served [Ne ii] fluxes and the star/disk properties 
such as infrared excess emission, stellar UV or X- 
ray flux, and mass accretion rates. 

First we explore any correlation between the 
continuum emission in the vicinity of the [Ne ii] 
line and the [Ne ii] line luminosities (see Fig. O. 
Although there is no obvious trend between the 
plotted quantities, it is interesting that the four 
detections cluster in a narrow range of [Ne ii] line 
luminosities (differences less than a factor of 2) 
and below 400 mJy of excess emission at 13 /im. 
The weaker mid-infrared continuum level (caused 
by an inner hole and/or grain growth) improves 
the line to continuum ratio in these relatively low 



^This estimate includes a factor of 2 uncertainty in the pre- 
dicted H I fiux, see Note to Table |4] 

''The coronal H II region lies on top of the slightly lower X— 
ray heated region, which is only partially ionized (~0.1 to 
1%). 
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Fig. 5. — Possible correlations of [Ne ii] line lu- 
minosities with disk properties. [Ne ii] detections 
are filled squares and right-most [Ne ii] limits are 
open squares with left pointing arrows. The num- 
bers correspond to the ID numbers of our targets 
as given in Table [T] Upper panel: Dereddened fiux 
between 12.85 and 13.15 ju m (as estimat ed from 



the low-resolution spectra. iBouwman et al.ll2007l ) 



times distance squared versus the [Ne ii] line lumi- 
nosity. To deredden the fluxes we used the visual 
extin ctions in Tablep] and the mean extinction law 
from iMathis I ( 1990l ). Lower panel: Excess emis- 
sion at 13 /xm versus the [Ne ii] line luminosity. 
The excess emission is computed from the differ- 
ence between the observed IRS low-resolution flux 
between 12.85 and 13.15 /jm and the best fit K u- 
rucz model atmospheres (see Meyer et al.l l200() for 
details on the procedure) in the same wavelength 
range. Excess emission has been dereddened as 
explained above. In both plots the errorbars in 
the y-axis are few percent of the plotted values 
and thus smaller than the used symbols. Low in- 
frared excess emission makes it more favorable to 
detect [Ne ii] emission lines. 
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Table 4 

Magnetospheric accretion model parameters and resulting mass accretion rates (M*). 



Source 








i 


log(M») 


Predicted Flux(H l) 




[Rq] 


[kms^i] 


[K] 


n 


[Meyr-i] 


[10"^^ ergs-^cm"^] 


RX Jllll. 7-7620 


2.0 


25 


12000 


75 


-9.3 





PDS 66 


1.0 





9000 


45 


-8.3 


2.4 


HD 143006 


1.5 


10 


12000 


10 


-9.3 


1.3 


[PZ99] J161411. 0-230536 


2.5 


45 


12000 


65 


-9.5 





RX J1842. 9-3532 


1.0 


25 


10000 


45 


-9.0 


0.48 


RX J1852. 3-3700 


1.0 


25 


12000 


35 


-9.3 


0.56 



Note. — Tmax is the maximum value of the adopted temperature distribution of the accretion column. 
All models are calculated for a solar mass star and for outer magnetospheric radii between 2.2 and 3 ii« . 
The last column gives the predicted flux in the H l(7— 6) transition at 12.37 fim. Uncertainties in the mass 
accretion rates are no more than a factor of 3—5. By varying by a factor of 5 the models give a range 
of H I fluxes within a factor of ~2. 




Fig. 3. — Observed and model Balmer line profiles (solid and dashed lines, respectively) for the sources 
RX J1852.3-3700, [PZ99] J161411.0-230536, and RX J1842.9-3532. Model parameters are listed in Table[ll 
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HP 143006 




Fig. 4. — Observed and model Balmer line profiles (solid and dashed lines, respectively) for the sources 
HD 143006, RX Jllll. 7-7620, and PDS 66. Model parameters are listed in TableSl 



spectral resolution Spitzer observations and thus 
makes it more favorable to detect [Ne ii] emission 
lines. 

We also search for correlations with the disk 
structure, in particular with the disk flaring. The 
continuum flux emitted from the surface layer of 
the disk is proportional to the angle at which stel- 
lar radiation impinges onto the disk. This so- 
called grazing angle becomes proportional to the 
disk flaring at radial distances fro m sol ar-type 
stars >0.4AU (Chiang & Goldrcicli[l993). Such 
distances are probed by dust emission at mid- 
infrared wavelengths. Thus, we can use the ratio 
of fluxes at two different mid-infrared wavelengths 
to trace changes in the disk flaring. We chose as 
reference flux that at 5.5 /xm, which is the short- 
est wavelength covered by the IRS low-resolution 
modules, and computed the flux ratios at 13, 24 
and 33fim (see Fig.|6]). Larger flaring is indicated 
by higher ratios of long wavelength continuum to 
short wavelength continuum flux from the dust 
disk. Source RX J1852. 3-3700 is pecuhar in having 
small excess at 13 /im and large excesses at 24 and 
33 /xm - The SED of this object (jSilverstone et al 



GM Aur and DM Tau (|Calvet et all 120051) . We 
do not see any obvious correlation between disk 
flaring and [Ne ii] line luminosities. 

Interestingly, we find trends with the stellar X- 
ray luminosities (Fig. [7]) and with the mass ac- 
cretion rates from the Balmer profiles (Fig. [5]). 
Stellar X-ray luminosities are computed from the 
ROSAT All-Sky Survey count rates and hardness 
ratios (see Tabled]). We did not attempt to correct 
for the interstellar extinction because the X-ray 
spectrum of our targets is not well-known from 
the ROSAT observations alone. Nevertheless, the 
largest difference in visual extinction (0.7 mag be- 
tween HD 143006 and RX J1852. 3-3700) trans- 
lates into less than 0.1 dex difference in log(ix), 
well within the estimated errors and variations due 
to stellar activity. This indicates that any ob- 
served trend in our sample can be only marginally 
influenced by the correction in interstellar extinc- 
tion. The best (minimum chi-square) flt to the 
[Ne ii] detections in Fig. Ogives the following cor- 
relation between the [Ne ii] and the X-ray lumi- 
nosity: 



20061 ) is reminiscent of transition systems such as 



-^[Nc II] oc 



1.2(±0.7) 



The slope of the relation is not constrained but 



11 



2.0 



„ 1.5 



6:^ 



1.0 

12 
10 



4 
2 
28 

15 

10 

ClT' 

5 



-•5- 



oL_ 

27.8 



«)2 
4 ., 1 



28.0 



28.2 28.4 
Log(LNeii) [erg/s] 



28.6 



Fig. 6. — Possible correlations of [Ne 11] line lu- 
minosities with the disk flaring. [Ne 11] detections 
are filled symbols and right-most [Ne 11] limits are 
open symbols with left pointing arrows. The disk 
flaring is represented by the flux ratio at 13 (upper 
panel), 24 (middle panel), and 33 /im (lower panel) 
over the reference flux at 5.5 /im, that is the short- 
est wavelength covered by the IRS low-resolution 
spectra. More flaring is indicated by higher ra- 
tios in the figure. There is no obvious correlation 
between disk flaring and [Ne 11] emission. 



it is positive within 1.5a. To further evaluate 
the likelihood of a positive correlation, we gen- 
erate 100 normally-distributed points around the 
four [Ne 11] detections and flt the data using six 
different linear re gression methods described in 
Isobe et al. (I1990I) . The analytical slopes as well 
the average slopes of 100 random resampling of 
the data using boot strap and jackknifc techniques 
(|Babu fc Feigelsonl[T99 2: Fcigclson & Babu 199j) 
^_ are all positive. The Pearson correlation coeffi- 

- cient is 0.24 meaning that we can be confident of 

- a positive correlation between L[Ne ii] and Lx at a 
I 95% confidence level. We perform a similar anal- 

- ysis to test the apparent anti-correlation between 
: the [Ne 11] luminosity and the mass accretion rate 

-1 in Fig. [H The ordinary least-squares fit to the 
four [Ne 11] detections results in: 



-^^[No II] oc Mi, 



-0.45(±0.36) 



'z The slope of the relation is negative only within 
; 1(7. The Monte Carlo approach in combination 
: with the six regression methods also indicates a 
- negative correlation at a 95% confidence level. 
: In summary, we have evidence for a positive 



correlation between L 



[Nc II] 



and Lx and a nega- 



- five correlation between L[Ne n] and M*. How- 
: ever we cannot quantify the slopes of these rela- 
2g ^ions because of the small number of detections 
and of the large errorbars in the measured and es- 
timated quantities. It is also important to point 
out that source 5 (RX J1842. 9-3532) has a large 
weight in determining these trends being the tar- 
get with the most different X-ray luminosity and 
mass accretion rate. Note that if the 2a detection 
from PDS 66 is real it would confirm the posi- 
tive trend between [Ne 11] luminosity and X-ray 
luminosity. A larger sample of sources spanning 
a wider range in X-ray luminosities and mass ac- 
cretion rates is necessary to constrain the slopes 
of the L[Ne ii]-ix and L[nc u]~M-i, relations. In 
the subsections below we discuss predictions from 
the X-ray and EUV models and compare them to 
the measured [Ne 11] fluxes and to the correlations 
presented in this Section. 

4.1. Plausible Disk Origin for the Ne+ 
emission 



Acc ording to the model proposed bv lClassgold. Najita &: Igeal 
POOTI ). detectable [Ne 11] emission can be pro- 
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Fig. 7. — Line luminosities (filled symbols) and upper limits (open symbols with downward arrows) for the 
[Ne ii] transition versus the s tar X-ray luminosity as given in Table [H The two open stars represent the two 
extreme thermal models from lGlassgold. Naiita fc Igeal (|2007l ): the upper model is when mechanical heating 
(accretion) dominates (Xray-M), the lower model is when X-ray heating dominates (Xray). The ordinary 
least-squares fit to the four detections suggests a positive correlation between L[nc ii] ^-nd Lx- 
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duced by the atmosphere o f a disk both ionized 
and h eated by stellar X-rays (jGlassgold. Najita fc Igea 
20041 ). Ne ions, primarily Ne"*" and Ne^"*", are gen- 



erated through X-ray ionization and destroyed by 
charge exchange with atomic hydrogen and radia- 
tive recombination. Because the X-rays emitted 
by young stars have a characteristic energy similar 
to the K-edges of Ne and Ne+ at ~ 0.9 keV, they 
are efficient in producing Ne ions. The warm disk 
surface region (T ^ 4000 K) generated by X-ray 
irradiation extends out to large radii (~ 20 AU). 
As a result, fine structure transitions of Ne+ 
and Ne^+, which have excitation temperatures 
of ~ 1000 K, can be produced over a vertical col- 
umn of warm gas of 10^^ — 10^" cm~^ and over 
a large range of disk radii. These circumstances 
lead to a significant flux of [Ne ii] at 12.81/im and 
[Ne III] at 15.55^m. 

For a typical T Ta uri disk dP'Alessio et"al 



119991) located at 140 pc, iGlassgold. Naiita fc Igea 
(|2007l ) predict that a stellar X-ray luminosity of 
log(Lx) = 30.30 erg s~^ generates a [Ne ii] flux of 
6.22 X 10~^^ ergs~^ cm~^ when accretion related 
processes are unimportant in heating the disk sur- 
face. The [Ne ill] flux is estimated to be approxi- 
mately 10 times lower. The X-ray luminosity as- 
sumed in the Glassgold et al. model is close to 
that of RX J1842.9-3532 (see Table [Hand Fig. [71) 
and the predicted [Ne ii] flux is similar (within 
a factor of 1.5) to the value we report in Table [3] 
Unfortunately, the upper limits on the [Ne ill] lines 
(see Table [3]) are not stringent enough to test the 
model predictions. 

In the model proposed by HoUenbach & Gorti 
(2007, in preparatio n), EUV photons fro m the 
stellar chromos phere ([Alexander et al. 2005 ) and/or 



stellar cnromos pnere l | Aiexander et ai.l zuuai ana/or 
from accretion ([Matsuvama et al .120031:1 iHerczeg et al 



20071 ) create an H ii-region like ionized layer on 



the surface of young circumstellar disks. The EUV 
photons incident on the disk are mostly absorbed 
near the so-called "gravitational radius" (~ 10 AU 
for a 1 M© star) where the local thermal speed of 
the 10^ K ionized hydrogen nuclei is equal to the 
escape s peed from the gravitational potential (e.g. 
Illollenbach et al.[[l994l) . The numerical model for 
the EUV-heated disk surfaces includes a calcu- 
lation of the diffuse field due to hydrogen and 
helium recombinations, and ionized gas chemistry 
comprising photoionizations, recombinations, and 
charge exchange reactions (HoUenbach & Gorti 



2007, in prep.). In the case of a soft EUV spec- 
trum (e.g. a blackbody at 4 x lO^'K), photon 
luminosities of ^ 10**^ erg s^^ produce a [Ne ii] 
line luminosity of ~ 10^^ L©, which corresponds 
to a flux of 2 X 10" erg s~^ cm~^ from a source 
at 140 pc. Harder EUV spectra can produce more 
doubly ionized Ne in the disk atmosphere which 
may result in stronger [Ne ill] lines in comparison 
to the case of soft EUV spectra (HoUenbach & 
Gorti 2007, in prep.). 

With the assumption that Ne atoms are ion- 
ized only by stellar EUV photons and that the 
sources have a soft EUV spectrum, one can use 
[Ne III lines as an indirect tool to estimate EUV 
fluxefl which are unconstrained for the majority 
of the stars. The [Ne ii] fluxes we report in Ta- 
ble[3lconvert to EUV fluxes between 2.6x10**^ (for 
RX J1842.9-3532) and 4.4xl0'*i (for RX J1852.3- 
3700) photons s~^. Such ionizing rates seem plau- 
sible for ~ 5 Myr old TTSs like our targets. For 
comparison [Herczeg et al.l (|2007() calculate ~ 10^^ 
ionizing phot ons s~^ for the ^10 Myr old TW Hya 



system while Alexan der et al.l (|2005l ) estimate a 
wide range of ionizing fluxes ~ lO''^ — lO^'' photon 
s"^ for a sampl e of five classical T Tauri stars. The 
estimates from [Alexander et al.l (|2005[ ) are based 
on modeling UV emission lines such as C iv and 
have only an order of magnitude accuracy due to 
model uncertainties and more significanly to un- 
certainties in the reddening. Nevertheless, if X- 
ra ys contribute to ionize Ne atom s, as proposed 
bv [Glassgold. Naiita fclgeal ^M), 

our EUV es- 
timates from [Ne ll] lines can be only taken as 
upper limits. 

In summary, both the X-ray and the EUV mod- 
els can reproduce the observed [Ne ii] line lumi- 
nosities with star/disk properties that are plausi- 
ble for our targets. Can the same models explain 
the H i(7-6) flux from RX J1852. 3-3700? HoUen- 
bach & Gorti estimate that Lh i ~ 10"'^L[no ii] if 
all the H i(7-6) emission originates from the ion- 
ized disk surface. They show that if X-rays dom- 
inate the [Ne ii] emission, then the H i luminosity 
should be even lower because the gas is mostly 
neutral. This demonstrates that neither the EUV 
nor the X-ray models can account for the observed 



'in the EUV model the [Ne ll] luminosity is directly propor- 
tional to the EUV luminosity of sources with soft EUV soft 
spectra 
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Fig. 8. — Line luminosities (filled symbols) and upper limits (open symbols with downward arrows) for the 
[Ne il] transition versus the stellar accretion rates given in Table HI We overplot a factor of 4 uncertainty 
in the mass accretion rates, representative of the uncertainties we estimate from the models (see Sect. 13.21) . 
The ordinary least-squares fit to the four detections suggests an anti-correlation between iv[No ii] ^-nd M^,. 
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H I flux, which is only a factor of ~2 lower than 
the [Ne ii] flux. We conclude that the H i(7-6) 
transition is not associated with either the X-ray 
excited gas or the EUV excited gas in the disk. 
Because magnetospheric accretion flows can ac- 
count for at most '^30% of the observed H i(7-6) 
flux (Sect. [321), accretion shocks and/or the stellar 
chromosphere are left as possible major sources of 
emission. 

4.2. Predicted Correlations and Observed 
Trends 

In this subsection we discuss whether the trends 
identified in Sect. 2] are consistent with predictions 
from the X-ray and EUV models. 

First we consider the correlation between 
the [Ne ll] luminosity and the X-ray luminos- 
ity. The X-ray model predicts that the ion- 
ization fraction of Ne atoms is proportional to 
the square root of the stellar X-ray luminosity: 
x{Ne+)lx{Ne) oc lI^'^ (eqs. 2.9 and 2.10 from 
Glassgold. Naiita fc Igeal l2007l ) . The integrated 



flux of the [Ne il] line further depends on the ex- 
citation of the line, which is a function of the disk 
temperature and electron density. Further model- 
ing by Meijerink, Glassgold & Najita in prep, indi- 
cate that there is a close to linear relation between 
-^[Ne II] and Lxj which is consistent with the posi- 
tive correlation suggested by our data. This trend 
could be also consistent with the EUV model. In 
the case of a soft EUV spectrum, most of the Ne 
in the disk surface is Ne"*" and the [Ne ii] lumi- 
nosity is predicted to be directly proportional to 
^Euv (HoUenbach & Gorti 2007 in prep.). What 
is the relation between Xeuv and LyJ Emission 
lines such as C IV and O vi tracing the EUV 
emission are found to scale with the X-ray flux 
with power la w indices close to 0.5 in ^Gyr ol d 
sun-like stars (jAvred [l997i iGuinan et"all l2003l ). 
If these relations hold for younger sun-like stars 
then the EUV model would also predict a positive 
co rrelation between Lj^e ii] and Lx- 

Glassgold. Najita k. Igeal ( 2007 ) also explored 



the effect of accretion heating on the [Ne ii] flux. 
Accretion at a level of 10~® yr~^, which is typ- 
ical for a classical TTS, results in a factor of 2 
higher [Ne ii] flux than that produced by a non- 
accreting TTS surrounded by the same circumstel- 
lar disk (see also Fig. [7]). This suggests a mild pos- 



itive correlation between i[Nc ii] and M^,, which is 
different from the one we see in our data (Sect. H 
and Fig. [8]). However a detailed comparison with 
the models needs to take into account other factors 
in addition to stellar accretion such as disk flaring 
and inclination. In the case of the EUV model, 
stellar winds greater than ^ 10^^" Mq yr~^ could 
drastically reduce the EUV photons reaching the 
disk (HoUenbach & Gorti in prep.). Given that 
the ratio of mass outflow rate to the mass ac- 
cretion rate is close to ~ 0.1 f or most classical 
TTSs with a large uncert ainty ( Shu et al. 2000l : 



White fc Hillenbrandl l2004[ l. these wind rates cor- 
respond to Mv, > IQ-^Moyr-i. Thus, the EUV 
model predicts significantly lower [Ne ii] luminosi- 
ties for classical TTSs simply due to their higher 
accretion (and thus wind) rates. 

Finally, disk flaring is also expected to play a 
major role in the X-ray model in that more flaring 
facilitates the penetration of stellar X-ray in th e 
disk atmosphere ( Glassgold, Najita fc Igeai2007l ). 
On the contrary, the [Ne il] emission produced by 
EUV photons is expected to be independent of the 
disk flaring (HoUenbach fc Gorti in prep.). Quanti- 
tative predictions from both models are necessary 
before attempting any comparison with the empir- 
ically derived signature of disk flaring presented in 
Fig. El 

5. Discussion 

We have shown in the previous Sections that 
both the X-ray and EUV models can account for 
the observed [Ne ii] fluxes suggesting that the 
emission originates in the hot disk atmosphere. 
Other sources of [Ne il] emission that we have not 
discussed so far might be a jet or an outflow. In 
this case, the jet/outflow needs to provide enough 
energy to ionize Ne atoms, i.e. about twice the 
energy to ionize atoms like S and N whose for- 
bidden ]iiiBS_aje_sometnn^ detected in classical 
TTSs (iHartiean et a l.l ll995h . Our sources however 
have very weak or absent mass loss signatures and 
accretion rates at least a factor of 10 lower than 
those of typical TTSs (see Sect. |321)- Thus, we 
regard the possibility of a jet/outflow dominating 
the [Ne ii] emission as very unlikely for our targets. 

It is possible to conflrm that [Ne ii] emission 
originates from the disk atmosphere by spectrally 
resolving the detected lines. High-resolution spec- 
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troscopy could also enable the detection of the 
weaker [Ne iii] line that can be used to further 
test the model predictions. Emission arising from 
disk radii out to ~ 10 AU around a solar mass star 
would produce a line with ~ lOkm/s width, about 
40 times smaller than the spectral resolution of 
the Spitzer IRS. Current ground-based spectro- 
graphs with spectral resolutions R > 30,000 at 
12/xm such as VISIR/VLT and TEXES/Gcmini 
should be able to resolve such emission lines. In 
addition similar observations of the H i(7-6) line at 
12.37 /im could confirm that the H i line originates 
in a different region from the [Ne ii] line. 

However high-resolution spectroscopy alone is 
not sufficient to identify whether X-rays or EUV 
photons dominate the ionization of Ne atoms. In 
fact both the X-ray and EUV model predict a 
similar extension of the [Ne ii] emitting region. 
Constraints for the EUV model can come from in- 
dependent estimates of the stellar EUV flux. Be- 
cause EUV emission is produced by gas at tem- 
peratures between 10^ — 10^ K, an estimate of the 
emitting plasma can be achieved by using gas lines 
such as the O vi tracing gas at intermediate tem- 
peratures. An important parameter in the X-ray 
model is the input X- ray luminosity and spec- 
tral te mperature which Glassgold. Naiita fc Igeal 
(j2007l ) ta ke from solar-like young stellar objects 
in Orion (jWolk et al.ll2005r ). It will be important 
to explore the effect of different X-ray luminosi- 
ties and spectra to better compare model predic- 
tions with observations. The anti-correlation be- 
tween [Ne ii] luminosity and mass accretion rate 
(if confirmed by a larger sample of sources) could 
be used to identify the dominant ionization mech- 
anism for Ne atoms. Detections of [Ne ii] emis- 
sion lines from classical TTSs, older TTSs, and 
young brown dwarfs will help sample the trend to 
mass accretion rates ranging from ~ 10^^ down 
to - 10"i°MQyr-i. 

None of the opt ically thin systems we published 



Pascucci et al.l f)2006.) exhibits [Ne ii] emission 



at 12.81 /Ltm. The majority of these sources have 
3(7 line flux upper limits in the H2 S(2) line (that 
are representative for the [Ne 11] lines) about an 
order of magnitude lower than the [Ne 11] fluxes 
we detec t in the four optically thick disks (see Ta- 
ble 4 in IPascucci et al. 20061 ). In addition, these 
optically thin systems span a wide range in X- 
ray luminosities that covers the X-ray luminosi- 



ties of the sources with detected [Ne 11] emis- 
sion. This indicates t hat th e optically thin sam- 
ple in Pascucci et al.l ( 2006h lack enough hot gas 
to emit detectable [Ne 11] emission. The EUV 
and X-ray models of optically thick dust disks 
as presented above can reproduce the observed 
[Ne 11] emission lines from only ^ 6 x 10^^ Mj 
of gas within 10 AU or ~ 3 x lO^'^Mj of gas 
within 20 AU, respectively. Thus, [Ne il] non- 
detections in optically thin systems might be used 
to set even more stringent gas mass upper lim- 



i ts th an t hose we reported in iHoUenbach et al 



( 2005 ) and Pascucci et al. ( 20061 ) using other mid- 
infrared gas transitions. However, a correct de- 
termination of the gas mass upper limits requires 
detailed modehng of the star/disk properties. We 
plan to explore this issue in a separate contribu- 
tion. 

6. Summary 

To summarize, the main conclusions of this pa- 
per are as follows: 

1. We detect [Ne 11] emission lines at 12.81 /im 
in four out of six optically thick dust disks 
observed as part of the PEPS Spitzer Legacy 
program. The systems with [Ne 11] emission 
are characterized by weaker mid-infrared 
continuum (possibly the result of an inner 
hole and/or grain growth) compared to those 
systems where we do not detect [Ne 11] lines. 

2. We also detect a H i(7-6) emission line from 
RX J1852. 3-3700. Magnetospheric accretion 
flows can account for at most ~30% of the 
observed flux. This line is not associated 
with the gas emitting the [Ne 11] lines. Accre- 
tion shocks and/or the stellar corona could 
contribute to most of the observed H i(7-6) 
emission. 

3. The [No 11] line luminosity correlates with 
the stellar X-ray luminosity. We find an 
anti-correlation between the [Ne 11] luminos- 
ity and the mass accretion rate. The slopes 
of these trends are not constrained by the 
current data. 

4. Emission from Ne"*" is very likely arising from 
the hot surface of the disk. Both stellar X- 
rays and EUV photons can sufficiently ionize 
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the disk surface to reproduce the observed 
hne fliixes. 
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